Approximate analytic solutions of the diatomic molecules in the 
Schrodinger equation with hyperbolical potentials 

Sameer M. Ikhdair^'Q and Ramazan Sever^'Q 

^Department of Electrical and Electronic Engineering, 
Near East University, Nicosia, Northern Cyprus, Turkey 
^Department of Physics, Middle East Technical University, 06531, Ankara, Turkey 

(Dated: September 8, 2009) 

Abstract 

The Schrodinger equation for the rotational-vibrational (ro-vibrational) motion of a diatomic 
molecule with empirical potential functions is solved approximately by means of the Nikiforov- 
Uvarov method. The approximate ro-vibratinal energy spectra and the corresponding normalized 
total wavefunctions are calculated in closed form and expressed in terms of the hypergeometric 
functions or Jacobi polynomials P'lf'''^\x), where /i > —1, v > —1 and x G [— 1,+1] . The s-waves 
analytic solution is obtained. The numerical energy eigenvalues for selected H2 and Ar2 molecules 
are also calculated and compared with the previous models and experiments. 
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I. INTRODUCTION 



The problems connected with the molecular sturucture provide interesting and instruc- 
tive applications of quantum mechanics, since molecules are considerably more complex in 
structure than atoms. Two distict problems arise in connection with molecular structure. 
The first is to obtain the electronic wave functions and potential energy functions of the 
nuclear coordinates. This problem can be solved analytically only in the simplest cases. 
The second is to obtain the solution of the nuclear motion equation. In solving the second 
problem, the construction of a suitable potential function of a diatomic molecule is very 
important. It has been found that the potential-energy function for the lowest electronic 
states of actual diatomic molecules can be expressed by the Morse potential [1] : 



which has three adjustable positive parameters a, D and r^.. At r = re, it has a minimum 
value at zero and approaches D exponentially for large r. If ^ is somewhat smaller than Tg, 
it becomes large (but not infinie) as r — > 0. This potential has been the subject of many 
studies since 1929 [1]. Also, it is an important potential in the field of molecular physics 
describing the interaction between two atoms and has attracted a great interest for some 
decades [1,2]. 

During the past years progress has been made in the field of diatomic molecules and 
extensive use of the potential functions have been introduced [3,4]. At present the Morse 
potential is still one of the potential functions used most in molecular physics and quantum 
chemistry [5] . However, the Morse potential has few asymptotic inaccuracies in the regions 
of small and large r. Obviously, unlike the true Coulombic interaction among two atoms 
which approaches infinity, the potential is finite when the distance of two atoms approaches 
zero. Another inaccuracy is the replacement of the Van der Waals term by an exponential 
which makes the value of the Morse potential smaller than that from the experiments in 
the region of large r. To avoid these inaccuracies, many works have been carried out in that 
direction to improve Morse potential [6]. 

In 1986, Schioberg [7] suggested hyperbolical (empirical) potential functions of the form: 



where D, a, 5, and a are four adjustable positive parameter with D — D^/ [5 — cr)^ {D^ is 



Vuir) = -D [1 - exp [-a{r - Tg)]]^ , 



(1) 




(2) 
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the spectroscopic dissociation energy) . It is similar to Morse potential having the minimum 
value at the point 



and approaches D exponentially for large r. The functional forms (2) are significant in the 
limits of < (7 < 5 for V^{r) and cr < (5 < for V-{r). V+(r) becomes more and more 
asymmetric if one lowers the potential constants 5 and cr but V_(r) becomes more and more 
symmetric with increasing 8 and a. One can change the potential constants of V±(r) in such 
a manner that the energy eigenvalues change into all known energy values which one obtains 
from rigorous quantum mechanical or semiclassical solutions. These are the energy values of 
the Morse, the Kratzer, the Coulomb, the harmonic oscillator and other potential functions. 
The Morse function represents the connection point between V-^.{r) and V-{r) [7]. Unlike 
the Morse potential (1), the empirical potential (EP) function V+{r) approaches infinity at 
the point r = 0. In the region of large r, it is closer to the experimental Rydberg-Klein- 
Rees (RKR) curve than the Morse potential for some diatomic molecules. This potential 
is not only a better description for the potential energy of a pure molecular vibrational (= 
radial) energies but it represents also perfectly intermolecular interactions and includes with 
the Kratzer energy eigenvalues the Rydberg terms of an electron in an atom. Besides, the 
semiclassically calculated term values for 'vibrations' in a plane there are also the quantum- 
mechanical calculated term values which include these radial terms and zero-point energies 
of a 'rotation' in space. 

Recently, Lu [8] solved approximately the Schrodinger equation of diatomic molecules 
with the EP functions using the hypergeometric series method. Furthermore, rigoorous 
solutions of the Schrodinger eqauation are also obtained with a similar method for zero 
total angular momentum. Since there are no exact analytic solutions for the EP functions 
(Z 7^ 0), some approximation [8] was used to obtain the solutions. This approximation was 
employed in solving the rotating Morse potential for any i-states [9,10]. The ro-vibratinal 
energy eigenvalues of the EP functions were determined with a semiclassical (SC) procedure 
(the Bohr-Sommerfeld quantization condition) and a quantum-mechanical (QM) method 
(the Schrodinger equation) [8]. 

The NU method [11] and other methods have also been used to solve the non-relativistic 
and relativistic wave equations [12-14]. The purpose of this work is to solve the radial 
Scrodinger equation with any orbital angular quantum number I for EP functions V±(r) 




(3) 
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using a generalized model of the NU method derived for the exponential-type potentials 
like EP. In addition, we apply the analytic solution to obtain the ro-vibratinal energy states 
for selected H2 and Ar2 diatomic molecules using the relevant potential parameters and 
spectroscpic constants given in Ref. [7]. 

The present work is organized as follows. In section 2, wc present a parametric general- 
ization of the NU method holds for any exponential-type potential. In section 3, we obtain 
the analytic NU bound state solution of the Schrodinger equation with the EP functions for 
any Z-states. In Section 4, we calculate the ro-vibratinal energy states for selected H2 and 
Ar2 diatomic molecules. Section 5 contains the relevant conclusions. 



II. NU METHOD 



The Nikiforov-Uvarov (NU) method is briefly outlined here and the details can be found 
in [11]. It is proposed to solve the second-order linear differential equation by reducing it to 
a generalized equation of hypergeometric-type of the following form: 

where the prime denotes the differentiation with respect to z, a{z) and (^{z) are polynomials, 
at most second-degree, and r(s) is a first-degree polynomial. In order to find a particular 
solution of Eq. (4), we decompose the wavefunction R{z) as follows: 

R{z)^(P{z)yr,{z), (5) 

which reduces Eq. (4) to a hypergeometric type equation 

a{z)y';{z) + T{z)y'^{z) + Xyn{z) = 0, (6) 

where (f){z) is defined as a logarithmic derivative 

cf>'{z)/cf>{z)^n{z)/a{z), (7) 

and the other part yn{z) is the hypergeometric-type function whose polynomial solution 
satisfies the Rodrigues relation: 

A 



where An is a normalizing factor and p{z) is the weight function satisfying the condition 

ia{z)p{z)r^r{z)p{z). (9) 
The function tt{z) and the eigenvalue A required in this method are defined as 



niz) = I [a'{z) - r{z)] ± ^J\ [a'{z) -t{z)]' -a{z) + ka{z). (10) 

and 

X = k + Tr'{z). (11) 

Hence, the determination of k is the essential point in the calculation of 'k{z), for which the 
discriminant of the square root in Eq. (10) is set to zero. Also, the eigenvalue equation 
defined in Eq. (11) takes the following new form 

X^^ -nr'iz) -^n{n-l)a"{z), n = 0,l,2,---. (12) 

where 

t{z)=t{z) + 2t^{z), (13) 

and it's derivative is negative (i.e., t'{z) < 0) which is the essential condition for any choice 
of proper bound-state solutions. Finally, comparing Eqs. (11) and (12), we obtain the 
energy eigenalues. 

In this regard, we can derive a parametric generalization of the NU method valid for any 

central and non-central cxponcntial-typc potential [12]. We begin by comparing comparing 
the following generalized hypergeometric-type equation 

R"iz) + In'Tl ^'^'^ + r n ^ (-^1^' + - ^(^) = 0' (14) 

Z[i.-C3Z) [2;(1-C32;)J 

with Eq. (4), we obtain 

t{z) = ci - C2Z, (15a) 
a{z) ^z{l- C3Z) (15b) 

^{z) = -iiz^ + - ^3- (15c) 
Further, substituting Eqs. (15a)-(15c) into Eq. (10), we find 

1 /9 

'k{z) = C4 + c^z ± [{cq - C3/C+ _) z^ + (c7 + /:;+_) z + , (16) 
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where 

C4 = ^ (1 - Ci) , C5 = ^ (C2 - 2C3) , C6 = C5 + C7 = 2C4C5 - ^2, Cs = C4 + ^3. (17) 

The discriminant under the square root sign in Eq. (16) must be set to zero and the resulting 
equation has to be solved for k, it yields 

= - (c7 + 2C3C8) ± 2^^, (18) 

where 

Cg = C3 (C7 + C3C8) + C6. (19) 

Substituting Eq. (18) into Eq. (16) and then solving, we obtain the physical choice of 
parameters: 

7t{z) = C4 + C5Z - [{^/c^ + cav^) z - ^/c^ , (20) 
A;_ = - (C7 + 2C3C8) - 2^/^. (21) 
In addition, Eqs. (13), (15) and (20) provide the parameter 

r{z) = 1 - (c2 - 2C5) z-2 [(v^ + C3v^) z-,/^, (22) 

whose derivative must be negative: 

r'{z) = -2c3 - 2 (^ + C3v^) < 0, (23) 

in accordance with essential requirement of the method [11]. Solving Eqs. (11) and (12), we 
obtain the parametric energy equation: 

(C2 - C3) n + csri^ - {2n + 1) cg + {2n + 1) (^ + c^^) + cr + 2c^Cs + 2^^ = 0, (24) 

for the exponential-type potential under study. Let us now turn to the calculations of the 
wavef unctions. The weight function p[z) can be calculated by means of Eq. (9) as 

p{z) = z'^^l - c^zY^\ (25) 

and consequently the first part of the wavefunctions throughout the Rodrigues relation (8): 

yn{z) = P^-iO'-")(l - 2c3^), Clo > -1, cn > -I, z e [0, 1/C3] , (26) 
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where 

2 

cio = ci + 2c4 + 2^ - 1 > -1, cii = 1 - ci - 2c4 + — > -1, (27) 

C3 

2C32;) are the Jacobi polynomials. Also, the second part of the wavefunctions 
can be found from (7) as 

cf>{z) = z^'^l - c^zY'\ C12 > 0, ci3 > 0, (28) 

where 

C12 = C4 + ^ > 0, ci3 = -C4 + — (i/c^ - C5) > 0. (29) 
Hence, the general wavefunction (5) has the general form 

R{z) = MnZ^^^l - C3^)""Pi"i°'"")(l - 2C3^), (30) 

where Mn is a normalization constant. 

III. BOUND STATES OF THE HYPERBOLICAL POTENTIALS 

The Schrodinger equation for diatomic molecules with hyperbolical potential functions 
takes the form: 

{ + % - ^±(r)]} V'n.^(r, ^, ^) = 0, (31) 

where u = . The wavefunctions in the above equation could be separated to the 

following form [15-18] 

^l^nimir:0,^) = ^Rni{r)YU9,ip), (32) 

where Yi^{9, (p) is a spherical harmonic with angular momentum quantum numbers I and m. 
The substitution of Eq. (32) into Eq. (31) gives the following radial reduced wave function 
Rni{r) satisfying 

§i-'-^+%lEn.-Vdr)]]K,(r)=0. (33) 

where ^^^^^ is the centrifugal potential with the boundary condition that Ri{r) vanishes near 

the points r = and r oo. Furthermore, if / is not too large, the case of the vibrations 
of small amplitude about the minimum, we can then use the approximate expansion of the 
centrifugal potential near the minimum point r = Tg as [19]: 

± exp(— 2Q;r) 



1(1 + 1) + , , , ±exp(-2ar) 

2 2 S ^0 + ^ITZZ 7 — n N + ^2 

\ 1 T exp(— ictrj 



1 + exp(— 2Q;r) 



(34) 



where 



Ao = 1- 



1 =F exp(— 2Q;re 



1 =F exp{—2are 



- (3 + 2are) 



Ai = ±2 [exp(2are) T 1] < 3 



^2 = [exp(2are) T 1] 



1 =F exp(— 2Q;re 
2are 

1 ^ exp(— 2are) 
2Q;rp 



- (3 + 2are) 



3 + 2Q;re — 



1 =F exp(— 2Q;re 
2 are 

1 ^ exp(— 2Q;re) 



(35a) 
, (35b) 
(35c) 



and higher order terms are neglected. In fact, Eq. (34) is the approximate expansion of 
the centrifugal potential l{l + l)r~^ and is valid for all r ~ rg, the minimum point of V±(r) 
since r is not singular there. However, the expansion is not valid near the singularity point 
r — 0. Overmore, it is a good approximation for small vibrations around the equilibrium 
separation r — r^- 

Substituting Eqs. (2) and (34) into Eq. (33) and introducing a new variable z — 
± exp(— 2Q;r), where z e [±1, 0] for V±(r), we obtain the hypergeometric-type wave equation: 

z{\. — z) z^[l — z) 



-[Kl + S^-Qi- 



(1 - zY 

+ {2Kl -Qi)z- K 



Rnl{z) = 0, 



(36) 



where Rni{z) = Rniir) and 



[0 - a) + , ^ ^ Ao 



s. 



2a^n? 
2j2D 



2a''h?' 



(37a) 
(37b) 

(37c) 



By comparing Eq. (36) with Eq. (14), we obtain specific values for the set of constant 
parameters given in Section 2: 

Cl = C2 = C3 = 1, C4 = 0, C5 = C6 = X^; + Sf - Ql, 

cr = -2Kli + Ql, cs = K^i, cg = Sf, 

Cio = 2X„;, Cii = 25";, Ci2 = K^l, C13 = 5"; + ^, 

= + sf - - 1 ^2 = 2X^, - Ql, ^3 = ^n^- (38) 



By using Eqs. (20)- (22), we find the following physical values: 



7r(^) = Kni + Kni + Sij z, (39) 

k^-Qi- 2KniSi, (40) 

r{z) = 1 + 2X„, - 2 (1 + Kni + Si) z, (41) 

where r'{z) = — 2 (1 + K^i + Si) < is the essential condition for bound-state (real) 
solutions. In addition, the energy equation can be found via Eq. (24) as 



Sf - Qi 



{Si + n + lY 



Si + n + ^ 



(42) 



The energy eigenvalues are obtained as follows. 



E^i ^D{5-af + 



l{l + l)h^ 



2/xr, 



2pL 



Sf-Qi-\-{Si + n + iy 



, n,/ = 0,l,2. 



Si + n+^ 

where n and I signify the usual vibrational and rotational quantum numbers, respectively. 
Thus, the ro-vibrational energy spectrum takes the following explicit form 

l{l + 1)^2 



De + 



2/xr, 



-Ao 



2/x 



&5 + S3(A.-Ao 



2^^? ^2 I i(i±i) 



K'+l) A I 1 



2m£^2 + 



(43) 



which is identical to Eq. (28) of Ref. [8]. Its important to mention that very similar 
expressions to the above expression for the energy states have been found over the past 
years for the hypebolical (exponential-type) potentials with 6 is being set equal to one in 
Eq. (2) (cf. [19-21]). Very recently, a new improved approximation [21,22] for the centrifugal 
potential term l{l + l)/r'^ was used different than the ones used in Ref. [8] and the one which 
is commonly used in literature by Ref. [23]. 

Let us now turn to the calculations of the corresponding wavefunctions for the EP func- 
tions. Thus, referring to the parametric generalization of the NU method in Section 2, the 
weight function in Eq. (25) takes the form 



p(z) = z'^^^'{l-z) 



\25, 



(44) 
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which gives the first part of the wavef unctions in Eq. (5) as 
Also, the second part of the wavefunctions (28) can be found as 



(45) 



(46) 



Hence, the unnormahzed wavefunctions (30) are being expressed in terms of the Jacobi 
polynomials as 

Rnliz) = KlZ^'-^il - z)^' + ipi2i.-2S.)(i _ 2,)^ (47) 

where M^i is a normahzing factor and p(^^ni,2Si) ^^-^^ _ ^ (^^"^+^^" ^Fi(-n, 2Kni + 2Si + 



n+ 1, 2Kni + 1; z) with (m), 



(m+n— 1)! 
(m-1)! 



is Pochhammer's symbol. 



Hence, the total wavefunctions of the EP are 

^ (r, e, (fi) = Afni- [±exp(-2Q;r)]-^"' [1 - ± exp(-2Q;r)]^'+5 
r 

X pi^Kr,i,2Si)(^^ _ ±2exp{-2ar))Yimi9,ip). 

where the normalization constants Afni are determined in Appendix A. 

Let us find the vibrational energy states for the s-waves {I — 0) from Eq. (43) as 



(48) 



2/1 



E„ — Dp — 



21 2 



, n = 0, 1,2, • • • ,n„ 



(49) 



where nmax is the number of bound states for the whole bound spectrum near the continuous 
zone, njnax is the largest integer which is less than or equal to the value of n that makes the 
right side of Eq. (49) vanish, that is, 

1 



-a6- 



(50) 



Thus, ninax cannot be infinite {En^^^^^ = Dg), which is reflected in the above condition. 
Furthermore, the corresponding wavefunctions for the s-waves can be easily found from 
Eqs. (37) and (48) as 

V'±(r,^,(^) =^fn^ [±exp(-2ar)]''" [1 - ± exp(-2Q;r)]^+^ 

X P^2^-2.)^^ _ ±2exp(-2ar))ro,o(^,¥'), (51) 

where the normalization constants jV„ are determined in Appendix A. 
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IV. APPLICATIONS TO DIATOMIC MOLECULES 



In this section, we calculate the energy states for the two selected H2 and Ar2 diatomic 
molecules using Eqs. (35) and (43). The spectroscopic constants of these two molecules are 
given in Table 1. The vibrating ground state energy eigenvalues (in cm"^) for the H2 
molecule in the EP functions V+{r) are found by means of parametric generalization version 
of the NU method for the potential parameters given in Table 2. Our numerical results 
obtained in the present NU model are listed together with the analogous numerical results 
obtained by using SC procedure and a QM method mentioned in Ref . [7] for various potential 
parameters. Obviously, as shown in Table 2, the results obtained in the present model are 
in high agreement with those obtained by QM, however, the SC procedure is proportionally 
different. Therefore, the differences between our results and SC procedure are less than 0.01 
cm'^, i.e., they are negligible because of these approximations: 1 a.m.u = 931.502 MeV/c^, 1 
cm-^ = 1.23985 x 10"^ eV and he = 1973.29 eV.A° [28]. The second application is applied to 
Ar2 molecule. We confine our study to calculate the ro- vibrating energy states for the V+(r) 
potential using the the following potential parameters: a — 25.23, 5 — 41.75 and a — 0.6604 
{A°)~^ [7] together with the parameters given in Table 1. The splittings of the energy states 
of s- waves E+ = E^{n ^ 0) — E^{n = 0) obtained by the NU method and SC procedures 
arc presented in Tabic 3. The present results AE^{NU) from NU method and AE^{SC) 
obtained from the SC procedures are also compared with four-different experimental results 
labeled by AE{a),AE{b), AE{c) and AE{d) taken from Ref. [7]. It is clear from Table 3 
that our results are very close with the expremimentally determined values as well as the SC 
procedure results. Finally, the approximated ro- vibrating energy states of the V+(r) given 
in Eq. (2) for the Ar2 and H2 molecules are also calculated for the I ^ case. Table 4 shows 
the energy levels for vibrational (n = 0,1,2,3,4,5) and rotational (/ = 0,1,2) quantum 
numbers. 



V. CONCLUSIONS 



In the present work, we have used a parametric generalization version of the NU method 
derived for any exponential-type potential to obtain the approximate solutions of the 
Schrodinger equation with any orbital angular momentum quantum number I for the hyper- 
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bolical (EP) functions. This method is systematical and efficient in finding the ro-vibrating 
energy states of a diatomic molecule and the normalized wave functions expressed in terms 
of the Jacobi polynomials. Obviously, the energy eigenvalues equation given in Eq. (43), 
considering the approximation in Eq. (35), is identical to those obtained by functional anal- 
ysis method [8]. The analytical result is tested in calculating the ro-vibrating energy states 
of the H2 and Ar2 molecules. Comparisons with the results of previous methods SC and 
QM for the ground state {n — I — 0) show that our calculations are in high agreement with 
those experimental results for H2 and Ar2 for low values of the ro-vibrating energy states, 
i.e., 1 — 0,1, 2. These systematical procedures could be useful for other molecular potentials 
in calculating their higher or lower approximated energy states [29-31]. 
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APPENDIX A: NORMALIZATION OF THE RADIAL WAVE FUNCTION 

In order to find the normalization constants J\fni, we start by writting the normalization 
condition: 

AA„f ^^J^ ^2i^--l(l - Zf^^+' ^p(2K„u2S0^-^ _ 2z)fdz. (Al) 

Unfortunately, there is no formula available to calculate this key integration. Neveretheless, 
we can find the explicit normalization constants Mni- For this purpose, it is not difficult to 
obtain the results of the above integral by using the following formulas [24-27,29-31] 
1 

y (1 - zf-^ z""-^ 2F1 (a, 7; az) dz = zF2 (i^, a, /9; + i^; T, «) , (A2) 



00 

and {a,b;c;z) = £(^±«±£) J. Following Ref. [24-27,29-31], we calculate 

p=0 

the normalization constants: 

^ [ r(2i^^, + l)r(25^ + 2) ^ (-!)'» (1 + n + 2{K^i + Si))^ T{n + m) 
[ 2ar{n) ^^m\{m + 2K^i)\T{m + 2{K^i + Si + l)) 

(A3) 
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where 

fni = 3F2 {'^Kni + m,-n,n + l + 2{Kni + Siy,m + 2{Kni + 5; + 1); 1 + 2Kni; 1) . 
Furthermore, the normahzation constants for the s-wave can be also found as 

^ \ t{2K + 1)T{2s + 2) ^ (-1)"^ (1 + n + 2{kn + s))^ T{n + m) 1 
2Q;r(n) ^ m\(m + 2kny.r(m + 2(kn + s + l)) 

m=0 

where 

fl'n = 3F2 {'^K + m, -n, n + 1 + 2(A;„ + s); m + 2(A;„ + s + 1); 1 + 2/c„; 1) , 

and 

where is given by Eq. (49). 
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TABLE I: The spectroscopic constants of the EP.for and Ar2 molecules [7]. 

Parameters H2 Ar2 
De (cm-1) 38281 99.55 
re (^°) 0.7414 3.759 
H (a.m.u) 0.50407 19.9812 



TABLE II: The EP parameters of the V+(r) and the ground state energy, (in cm ^) of the H2 

molecule. 

a 5 a{A°)-^ E+{SC) E+{QM) Present 

426.826 463.102 0.9327 2167.68 2168.93 2168.68 

47.294 102.341 0.6146 2153.69 2164.83 2164.45 

28.685 117.121 0.3826 2139.57 2157.69 2157.53 

21.250 213.212 0.1762 2124.29 2148.40 2147.53 
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TABLE III: Comparisons of experimentally calculated s-states energy transition values 
AEnfl{cm^^) for n / ^ n = together with the results of the SC procedure and the present 
NU method for the Ar2 molecule. 



n Present AE{a) AE{b) AE{c) AE{d) AE+{SC) 



1 25.808 


25.74 


25.49 


25.21 


25.56 


25.75 


2 46.079 


46.15 


45.63 


45.02 


46.00 


46.01 


3 61.472 


61.75 


60.70 


60.04 


61.32 


61.42 


4 72.536 


72.66 


71.33 


70.92 


71.52 


72.52 


5 79.733 


79.44 








79.79 


6 83.453 










83.59 


7 84.026 
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TABLE IV: Energy levels En,i{cm ^) for Ar2 and H2 molecules in V+(r) using the NU method. 
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